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Abstract 

Band structure of InSb thin films with < 100 > surface orientation is calculated using empiri- 
cal pseudopotential method (EPM) to evaluate the performance of nanoscale devices using InSb 
substrate. Contrary to the predictions by simple effective mass approximation methods (EMA), 
our calculation reveals that T valley is still the lowest lying conduction valley. Based on EPM 
calculations, we obtained the important electronic structure and transport parameters, such as 
effective mass and valley energy minimum, of InSb thin film as a function of film thickness. Our 
calculations reveal that the 'effective mass' of F valley electrons increases with the scaling down 
of the film thickness. We also provide an assessment of nanoscale InSb thin film devices using 
Non-Equilibrium Green's Function under the effective mass framework in the ballistic regime. 

PACS numbers: 73.40.Gk, 73.40.Rw, 75.70.Cn 
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I. INTRODUCTION 



Ultra-thin body (UTB) metal-oxide-semiconductor-field-effect-transistor (MOSFET) 
structure is a promising candidate for scaling MOSFET devices into the nanometer regime 



because it has the excellent attribute o 
the downscaling of device gate length 



suppressing various short channel effects caused by 
l|. Recently, there are also experimental and theo- 

2|-y] or new channel 



retical efforts to evaluate use of non-conventional channel orientation 
materials (such as Ge, III-V compound semiconductors, due to their small F-valley electron 
masses) [4-^ to improve the MOSFET performance. 

In particular, InSb is being considered as a new channel material {gI for state-of-the-art 
MOSFET devices because of its high bulk mobility (in fact, it has the highest electron 
and hole mobility among common 111-V semiconductors [3, 1^)- Theoretical investigation of 
InSb MOSFET devices based on a simple effective mass approximation (EMA) had been 
conducted in Q to predict its device performance limits . It was predicted that the F-valley 
energy rises rapidly under body quantization and the inversion charges in the thin film are 
transferred to the L-valley. Therefore the advantage of the high injection velocity from the 

n 

F-valley electrons is lost. Work in [9.] used bulk effective mass even in thin film regime. 
However, EMA may not be a reliable method in describing the electronic band structure 
of thin films due to unaccounted effects like band coupling and non-parabolic dispersion. 
We also found that the bulk effective mass can only describe the energy dispersion in a 
very small range of k space, which brings to question how accurately EMA can capture the 
quantization effect in thin film (see the discussions on Fig. 7 in this paper). 

In this work, we describe a more physically reliable and accurate model, the local empirical 



pseudopotential method (EPM) 10|, |n|, to study InSb thin films. Our results show that 



the charge transfer from the F to the L valley does not occur and the effective masses 
in transverse directions become larger than those in bulk material, which have a direct 
impact on the device transport properties. This paper serves to communicate in detail the 
models employed in our published conference paper 3|. In section II, we give a detailed 
description of the theoretical background of EPM. We will discuss introduction of a model 
potential to describe the atomic pseudopotential, which allows us to extend EPM to thin 
film calculations. In subsection II. B, we derive the spin-orbit coupling contribution to the 
matrix element of the Hamiltonian used in EPM calculation. In order to achieve accurate 
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results in EPM calculations of thin film band structures, we introduce two parameters to 
take account into the volume renormalization and the spin-orbit coupling. In subsection II. 
C, we discuss the methodology adopted to passivate the surface dangling bonds by using 
Hydrogen (H) atom bonding. In section III, we discuss the important features of InSb thin 
film electronic structure. We also obtained important electronic parameters, such as effective 
mass and valley energy minimum, of InSb thin film as a function of film thickness. In section 
IV, we provide an assessment of nanoscale InSb thin film double-gated MOSFET devices 
using non-equilibrium Green's Function under the effective mass framework in the ballistic 
regime derived from realistic thin film band structure reported in this paper. Finally, we 
give conclusions in section V. 



II. THEORETICAL BACKGROUND 




FIG. 1: (Color online) Supercell of InSb used for our calculation. It consists of a thin film layer 
and a vacuum layer in the quantization direction. 9 is the angle of the H-In bond as illustrated. 
The positions of In and Sb atoms are the same as in the bulk material, z direction is in < 100 >. 



The EPM has already been successfully applied in band structure calculations of metal. 



semiconductor or other materials 



ulate band 

structures and electronic properties of Si quantum dots 16|, ll7|], quantum wires 18|, Il9| . 
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and quantum films 20|, |2l| . In view of pedagogical clarity, we shall also give an overview 



of the essential theoretical aspects of EPM. Before going into the rigorous mathematical 
details, we shall first give an overall description of the method used in this work. The key 
aspects of applying the EPM to a thin film calculation are as follows; (1) a supercell (Fig. 
1) has to be specified prior to any calculation. A sufficiently large vacuum layer has to 
be explicitly included in the supercell so as to eliminate any wavefunction overlap between 
different thin film layers. (2) We have to determine the suitable atomic pseudopotentials 
for the Sb and In atoms independently for our thin film calculations. (3) The effect of 
spin-orbit coupling needs to be included because it is sufficiently large to affect the band 
structure significantly (see Table III) in InSb and many other compound semiconductors. 
(4) Hydrogen passivation of the surface dangling bonds is incorporated into model potential. 
This is very crucial treatment as the surface dangling bonds will lead to surface states in the 
band gap which will overwhelm the band structure of the thin film, especially in small band 
gap materials like InSb. Therefore we would also need to determine the pseudopotential of 
Hydrogen atoms. 



A. Pseudopotential Method 

1. Empirical Pseudopotential of Bulk InSb 

It is well-known that the band structure can be derived by solving a secular equation j^]. 
The matrix element of the pseudopotential (PS) is 

{PW, G|y^^|Piy, G') = J2 exp(-^(G - G') • r)UP{G - G'), (1) 

r 

where 

0'{G) = i / UPiv) exp(-zG ■ r)dr, (2) 

Q is the volume of the unit cell, is the number of unit cells, G and G' are reciprocal 
lattice vectors, Uf^lG) is the form factor, and \PW, G') is the plane wave basis function. 
In a single crystal, there are usually several basis atoms associated with each lattice point 
fj located at a position in a unit cell. These atoms have relative positions r from the lattice 
point at fj. The summation over r will yield us the atomic configuration information of the 
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unit cell. Note that, the bold letters and the letters with an overhead arrow denote vectors 
throughout this paper. 



Form factor (Ry) 


q^ = 


q' = 3 


g2 = 4 


g2 = 8 


g2 = 11 


q^ = 12 


ESAFF 




-0.858 



-0.20 
0.035 




0.032 


0.018 



0.034 
0.011 




0.013 


EMP 




-0.816 



-0.202 
0.0353 




0.0312 


0.0179 



0.03416 
0.01221 




0.0114 



TABLE I: Symmetry U^{q^) and anti-symmetry form factors U^{q^) in our calculations. Unit of 
q^ is (27r/a)^, where a is crystal lattice constant. For InSb, we use a = 6.47877A. Form factors 
derived from empirical model potential (EMP) are also given to establish EMP validity. 



Before we proceed with the calculation of InSb thin film band structure, we need to 
derive the atomic pseudopotential (PS) of In and Sb atoms respectively, or more accurately 
speaking, to construct an empirical PS. We can obtain important information about the 
atomic PS from bulk InSb band structure, since its calculated band structure can reliably 
be calibrated against experimental data. For bulk InSb, there are two atoms in a unit cell. 
Conventionally, bulk empirical PS (at special values, they are also known as Empirical 
Symmetry and Anti-symmetry Form Factors (ESAFF), as shown in Table I) can be defined as 
the summation and the difference of the individual atomic PSs of Sb and In. The symmetry 
and anti-symmetry structure factors are given respectively, as = cos{9') and = sin(6''), 
where 6' = AG = G — G'. It should be noted that if we define the symmetry and anti- 
symmetry form factors for InSb as = (1/2) (f/|;^ + Uf^^) and = (1/2) (f/|;^ - Uf^^), 
then the factor 1/Q appearing in Eq.(l2]) will be 1/flatomic, where flatomic is the atomic volume 
in a unit cell. 

By performing EPM calculation of the bulk band structure iteratively, a suitable set of 
!5SAFF is obtained such that it yields the 'correct' energy gaps at pertinent symmetry points 
71, Is] . This set of ESAFF is shown in Table I. Because of the spherical symmetry approxi- 
mation of the local pseudopotential, the first few lowest energy shells are at normalized 
equal to 3, 4, 8, 11 and 12, where q is wave vector in unit 271 /a and a is crystal lattice 
constant. Based on this set of ESAFF, we can devise a suitable interpolation scheme that 
allows us to reasonably predict a pseudopotential value at other values. This will then 



allow us to extend this EPM method to calculation of thin film band structure of arbitrary 
thickness. An essential assumption we are invoking is that the atomic pseudopotential is 
transferable and thus will not be changed from the bulk to a thin film. This assumption is 
reasonable because ab initio calculation has demonstrated that the bulk potential will just 
be changed in only one atomic layer at the interface for GaAs/AlAs superlattice 22|] and 
will be exactly the same for all other atoms. Hence, the potential will be bulk-like potential 
except at one atomic layer at the surface of the thin film. Therefore, the bulk atomic po- 
tential can be used in thin film calculation. However, further investigation reveals that we 
must add the renormalization factor to the pseudopotential to give the correct results for 
thin film as shown in the subsection II. A. 2. 
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FIG. 2: (Color online) (a) Three-dimensional Brillouin zone of InSb and illustration of the impor- 
tant symmetry lines and points F, A, L, X and U [7]; (b) The two-dimensional Brillouin zone 
obtained by projecting onto the xy plane. The densest shaded region is the first Brillouin zone. 



2. Thin Film and Empirical Model Potential (EMP) 

For < 100 > surface InSb thin film band structure calculation, we use a supercell as 
shown in Fig. 1. Fig. (2a) depicts the first Brillouin zone of bulk InSb and its subsequent 
projection onto the xy plane, i.e. the first Brillouin zone of InSb thin film in 2D, is shown 
in Fig. (2b). In this paper, band structure of InSb thin film will be plotted using the same 
coordinate system for 2D reciprocal vector space as illustrated in Fig. (2b). 

In bulk InSb, band structure is determined by the pseudopotential at several values 
as stated in previous section and the shells with < 3 are not of much importance. For 
thin film InSb, the reciprocal primitive vector will be shorter than that of bulk and thus the 
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radius of the lowest energy shell is smaller. This is the consequence of a larger primitive 
vector in real space in z direction for < 100 > film (see Fig. 1). Therefore, the shells with 
< 3 are also very important for band structure calculation and their effects must be 
properly predicted and included in our calculation. To account for them, we shall employ 
an empirical model potential that can give the continuous pseudopotential as a function of 



q. We use a model potential 



17 
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2l| 



PS 
atom 



as exp (a4g^j — 1 

where the subscript " atom" is used to distinguish the different atomic species' pseudopoten- 
tial. ai, 02, 03 and 04 are fitting parameters such that they yield the correct form factor as 
determined by ESAFF. 




q'[(2jt/a)'] 



FIG. 3: (Color online) Atomic model potential of Sb and In. X-axis and y-axis are in units of 
and Ry (Rydberg), respectively. 



'27r\2 
> a ' 



Parameters Sb In 

ai 0.2588 719470 

02 1.5832 2.0811 

as 1.9689 3813600 

04 0.7159 0.9116 

TABLE II: Derived parameters for model potentials of atoms Sb and In. 



To obtain a unique atomic model potential of In and Sb, we need eight constraints for the 
eight fitting parameters. Seven constraints come from the seven nonzero ESAFF parameters 
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at equal to 3, 4, 8, 11 and 12 shown in Table I. The remaining one constraint is imposed 



by InSb work function (WF), which gives one nonzero ESAFF at equal to 
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last constraint accounts for the surface property, which is important for a thin film 



This 



2l|. In 



fitting these parameters, it should be noted that the physical unit of in Eq. dH]) is (27r/a) 



and is in rydberg atomic unit. The best fitting parameters for the model potential are 
shown in Table II and its variation with q^ is shown in Fig. 3. From this figure, we observe 
that the model potentials of Sb and In atoms are in a reasonable range. Also, the fits to 
and U"^ are quite good as is clear from Table I. 



Energy (eV) 


Experiment 


Calculation 1 (ESAFF) 


Calculation 2 (EMP) 


Et 


0.17 


0.172 


0.168 


El 


0.68 


0.685 


0.704 


Ex 


1.0 


0.995 


1.034 


Eso 


0.8 


0.801 


0.802 


WF 


-4.76 


-4.760 


-4.223 



TABLE III: Experimental data [7', 'S\ of various energy conduction valley minima {Er, El and 
Ex)-, spin-orbit coupling {Eso) ^-nd work function (WF) are compared with theoretical result from 
empirical pseudopotential method. Calculation 1 and 2 employs the Empirical Symmetry and 
Anti-symmetry Form Factors (ESAFF) and Empirical Model Potential (EMP) respectively shown 
in Table I. Although the error of WF derived from EMP is large, the relative error is still tolerable. 



We perform EPM calculation of bulk InSb using this new model potential and reproduce 
the required energy band minima derived from experiment, as tabulated in Table III. We 
shall point out that the model potential of one of the species of an atom in a semiconductor 
compound is actually atomic model potential, i.e. the volume of unit cell is divided by 
the number of atoms in this unit cell. So for the thin film calculations, there will be a 
renormalization with a different unit cell volume under the assumption that the overall 
pseudopotential of atoms is the superposition of the local pseudopotential of all the atoms 
in the thin film. 
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I 
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FIG. 4: (Color online) Energy bands of Si < 100 > thin film (8 atomic layers Si film + 8 atomic 
layers vacuum with H passivation) calculated with ab initio method (using CASTEP, Local Density 
Approximation) and via empirical pseudopotential method with a model potential. Excellent 
agreement for the various conduction valleys minima is obtained. 

B. Hydrogen Passivation 

In thin film calculation, we must ensure that surface dangling bonds are correctly passi- 
vated with Hydrogen (H) so that there will be no surface states within the band gap . To 
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2l| and 



determine model potential for Hydrogen, we tried model potential of H in 
calculated the band structure of 8 atomic layers (atm) Si < 100 > thin film and compared 
the EPM results with those derived from ah initio calculation [9] (Fig. 4). The band 
structure obtained from both methods are very consistent, thus establishing the validity 
and reliability of this H model potential. In Si < 100 > thin film calculation, the renormal- 
ization of supercell volume is set as a ratio of the number of H atoms in the supercell and 
the number of Si atoms in the same unit cell. Also we assumed that the atomic volume of 
H is the same as that of Si. 



For InSb thin film, the bond length of In-H is set as 1.75A (from experiment in 24]) and 



Sb-H as 1.711 A 
1.67A in Ref. 



rom ah initio energy minimization calculation [23], which is very close to 
25(1). The angle between the bond of In-H and the top surface {0) (Fig. 1), 
should be determined carefully. Fig. (5a) shows the calculated direct InSb band gap (i?r) 
as a function of 9 for varying number of InSb atomic layers. In principle, the energetically 
stable 9 for the thin film system is when Ey is stationary with respect to 9. Er as a function 
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FIG. 5: (Color online) (a) The calculated InSb direct energy gap {Ey) as a function of H-In bond 
angle 9, plotted for InSb thin films of varying number of atomic layers (atm). The graph illustrates 
that there is a common 9 where the energy function is a stationary point, (b) The direct energy 
gap {Ey) at bond angle 9 where the energy function is a stationary point, plotted as a function 
of film thickness. An expression which fits these data points is given. It gives correct asymptotic 
behavior when film thickness tends to bulk, yielding bulk band gap value of 0.43 eV. 



of 9 shown in 
(similar to 



big. 
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'5a) confirms 9 = 0.77rad as a stationary value for 9 and a fitting formula 



2l|) affirms that Ey approaches the bulk band gap of 0.43 eV as the film 



thickness tends to the bulk limit as shown in Fig. (5b) (without spin-orbit contribution). 

For InSb thin film, we need to renormalize Hydrogen atomic potential derived from Si 
thin film calculation, V^'^(InSb) = (i^si/^inSb)^jy "^(Si) using relative volume of H atoms 
compared to the volume of InSb, where flsi and fiinsb are supercell volumes of Si and 
InSb respectively. Since the InSb thin film has the same supercell structure as Si thin 
film, its supercell renomalization is performed using their respective lattice constants as 

(fisi/^InSb) 



(asi/ CllnSb)^- 



C. Spin-Orbit Coupling (SOC) 

SOC treatment is incorporated by calculating the matrix element of the spin-orbit cou- 
pling term in the Hamiltonian expressed in a plane wave representation. It can be reduced to 
an expression similar as the pseudopotential form factor, known as the spin orbit form factors 



(SOFFs) 26l-l29| . In this section, we will clarify the formulation used for SOC treatment. A 



more mathematically detailed derivation shall be presented in the Appendix. 
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If we consider a supercell composed of different species of basis atoms at each lattice site, 
the spin-orbit Hamiltonian can be specified as 



H 



so ^ 

k's' ,ks 



■ {k' X fc)^5j(A k)\ 



J' 



(4) 



where is a parameter for j-th species, m is the number of species, A k = k' — k is the 
change in the wave vector, a^/g are matrix elements of Pauli matrices, s and s' denote the spin 
states; up spin or down spin. Structure factor is given as Sj{A k) = n~^^ exp{i A k-Tn-), 
where the summation is over nnum atoms of j-th species (see Appendix for more detail). 
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FIG. 6: (Color online) (a) Variation of B function with wave vector for Sb atom and In atom for 



n = 4: and I = 1. (b) Its relative difference P with wave vector, where P={B^j'-B^^)/{B^j'+B^^^ 



The parameters are determined by the B function defined in the appendix (see Eq. 
( 1A19I) ). The spin-orbit parameter t] (see Eq. ( 1A15|) ) is a property of isolated atoms, so 
it is independent of the crystal wave vector. The defined B function is important in the 
calculation of Xj and its behavior is depicted in Fig. 6. 

In Fig. 6, we plot the variation of the B function (i.e. Eq. flAlOp ) as a function of wave 
vector for Sb and In atoms in (a) and the difference of the B function between these two 
types of atoms in (b). For our plot in Fig. 6, we just show the outermost core state, i.e. 
n = 4 and I = 1 {p orbit). Because of the normalization, the B function tends to 1 as — )■ 0. 

Although the explicit form of the spin-orbit coupling matrix element is presented in the 
Appendix, there are some additional steps that must be taken into consideration in order to 
obtain the correct result. Thus, we shall briefly elaborate the procedures involved in SOC 
treatment. First off all, the SOFFs are obtained via a similar methodology as ESAFF; by 
iterating our EPM calculations until the correct energy gaps are obtained for bulk InSb. 
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FIG. 7: (Color online) (a) Energy band structure of bulk InSb calculated using empirical pseu- 
dopotential method. We have calculated the various effective masses for the different conduction 
valleys. The longitudinal mass (derived from the empty circles in (a) at L, X and F) for L and X 
valleys are found to be = 2.45mo and mx = 3.90mo, respectively. Lastly, the isotropic effective 
mass at F valley is found to be m-p = 0.016mo. (b) Effective mass of F valley (mr) calculated from 
the second derivative of energy dispersion. 

SOFFs have to be normalized with respect to supercell volume when used for a thin film 
calculation. Besides the usual volume renormalization, one has to account for another 
renormalization. Conventionally, bulk band structure of InSb is derived by considering 
a 2-atom unit cell in which the primitive vectors are non-orthogonal and the mutual angles 
between them are all 7r/3. In thin film calculation, a larger supercell is constructed in which 
the primitive vectors are mutually orthogonal. Hence this angular effect renormalization 
and volume renormalization must be considered in our treatment. 

The spin orbit Hamiltonian is divided into two parts; one part dealing with the angular 
part and the other the magnitude part. When the supercell is changed from bulk 2-atom 
unit cell to < 100 > thin film supercell, the volume renormalization will appear in the 
second part like what has been done on form factors of pseudopotential of Sb and In. For 
the angular part renormalization, we treat it as an unknown parameter which we derived 
using the following procedure. We removed the vacuum slab from the thin film supercell. 
This new supercell is therefore a unit cell for a bulk crystal. We then adjust the angular 
part renormalization parameter so that it reproduces the correct bulk band structure and 
spin orbit splitting (as calculated from a conventional bulk crystal unit cell). This procedure 
has to be repeated for each film thickness. After this remormalization, the vacuum slabs are 
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replaced again and the thin film band structure can be derived from EPM. 




III. FEATURES OF InSb ELECTRONIC STRUCTURES 

Fig. (7a) shows the calculated band structure of bulk InSb from our EPM. We derived the 
following effective masses: rriL = 2.45mo (the longitudinal mass of L valley), mx = 3.90mo 
(the longitudinal mass of A valley) and m-p = 0.016mo (the isotropic mass of F valley). 
The open circles in Fig (7a) are data points fitted to the bands minima by using a parabolic 
dispersion with effective masses as stated. Fig. (7b) shows the 'effective mass' in the vicinity 
of F valley (mr) calculated by taking the second derivative of energy with respect to the 
wave vector k. We note that the parabolic assumption for the energy dispersion at F is only 
valid for a very small k range. Consider a 1.3 nm InSb film, the wave vector spread according 
to uncertainty principle is ~ 0.04(27r/a). From the plot, the parabolic assumption for EMA 
is only applicable up to a A; range of ~ 0.01(27r/a). Hence, in the ultra thin film regime, 
a parabolic EMA is not a reasonable assumption and becomes a highly unreliable method 
for calculation of size quantization effect. In addition, the small band gap also entails a 
considerable amount of coupling between the conduction and valence bands, which render 
the uncoupled EMA approach unreliable. 

Fig. 8 shows the calculated band structure of thin film InSb from our EPM. These thin 
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FIG. 9: (Color online) (a) Variation of band gap with the number of atomic layers. (Top scale 
shows the thickness of thin film in nm in (a) and (b)). (b) The InSb thin film valence band 
maximum (VBM) and conduction band minimum (CBM) as a function of the number of atomic 
layers. Bulk values are approached at thicker films. 

n 

films still retain the direct band gap properties in contrast to EMA calculation [9] , which has 
L valley as the lowest lying conduction valley when film thickness is below 5 nm. Fig. (9a) 
shows variation of the energy minima at F, L and A valleys as function of film thickness. 
We also plot the fitting curves for the energy minima as a function of thickness x in this 
figure. The fitting formulas are 0.16802 + 150/(x + 'd.lf -^^, 1.03395 + 92/(x + h.hf-^^, and 
0.70439 + 94/(x + 5. 2)^'^^, respectively. There was no crossing of the energy minima down 
to 6 atomic layers, in contrast to predictions by EMA methods. In addition, we note that 
the lowest lying F valley is separated from the other valleys with a gap more than 0.3 eV, 
signifying that electrons are dominantly occupying F valley. Another starking contrast with 
results of EMA is the band gap. Although, the band gap is enhanced by quantization effect, 
this increase with reducing thickness is much slower than the well-known dr"^ behavior 
predicted by EMA in an infinite deep well model (For e.g. at F valley, see Fig. (9a)). 
However, EPM also yields the result that the conduction band quantization effect is larger 
than that of valence band in agreement with prediction by EMA (Fig. (9b)). From a logic 
device point-of-view, one would desire to have a larger band gap to curb the increasing 
band-to-band tunneling current with each new generation CMOS devices. The fact that 
InSb is a direct band gap material may aggravate the problem. However, if one can achieve 
a sufficiently thin film with ~ 8atm, which offers a band gap more than 1 eV, band-to-band 
tunneling current should still be tolerable. 
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FIG. 10: Energy contour plot for 8 atomic layers InSb thin film for the first conduction and first 
valence band are plotted in (a) and (b), respectively. Values of the energy are indicated on each 
contour line. Important symmetry points are also indicated. It is apparent that the T conduction 
valley is relatively isotropic in nature even in the ultra-thin film regime. 

Fig. 10 shows the contour plot of the conduction and valence band energy dispersion for 
8 atomic layers InSb thin film. We confirm that the F valley in thin film still retains its 
isotropy unlike the case of Si ^. Fig. 11 shows the Si and InSb thin film effective mass 
fitted from their 2D energy dispersion. InSb F valley is isotropic with increasing effective 
mass as film thickness is scaled down. Si F valley in thin film actually is originated from bulk 
A valley, projected onto the 2D k-space. We observe that anisotropy in Si becomes more 
prominent with decreasing of film thickness, with the TK direction effective mass diverging 
for each of the two degenerate bands. The larger TK effective mass in Si is for the lower 
band of the two degenerate bands, and the smaller TK effective mass corresponds to the 
top band. 

We have attempted to fit a parabolic dispersion to match the calculated energy dispersion 
to obtain the effective mass. In fact, we expect the effective mass to also increase with wave- 
vector as depicted in Fig. (7b). However, even if effective mass is not a rigorously derived 
quantity due to non-parabolic valleys, it serves as a very useful 'figure of merit'. Interestingly, 
the isotropic effective mass of InSb F valley increases with the decrease of film thickness. 
At ~ 1 nm of InSb film, the effective mass is ~ O.lmo- Increase in mr would retard the 
quantization effect in InSb thin film. The isotropy of the electron mass, which translates to 
isotropy of the electron transport property, should be advantageous as it affords engineer 
with more flexibility in orientating the n and p-MOSFET devices to yield the most optimum 
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transport direction on the same substrate. 
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FIG. 11: (Color online) Electron effective mass of Si and InSb F conduction valley. The inset shows 
the constant energy ellipsoid for Si projected onto the 2D k-space, with the F valley being doubly 
degenerate. 
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FIG. 12: (Color online) Ballistic current calculated for InSb and Si double-gated devices using 
simple effective mass approximation. 8atm and 24atm InSb/Si thin film devices are considered. 
Effective masses used for InSb are derived from a parabolic fit of the band structure from EPM. 
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IV. BALLISTIC LIMIT OF InSb NANOSCALE DEVICES 



We calculated the InSb thin film devices' ballistic current limit using the Non-Equilibrium 
Green's Function (NEGF) method 30|], under the framework of simple effective mass 
method. Since the energy dispersion for < 100 >-surface InSb thin film is relatively isotropic 
(Fig. 10), one can employ a decoupled 2D treatment to the problem and calculation is done 
in the framework of mode-space NEGF approach 30|. The effective masses used for the 
various InSb devices are as derived in Fig. 11. For Si devices, due to the anisotropy, we 
simply assumed an effective mass of 0.20mo. A double-gated device structure shown in Fig. 
12 is employed. Channel length of 20 nm is used. An EOT of 1 nm and metal gate is 
employed. No wave function penetration into oxide is assumed. N-type Source/drain and 
p-type channel are doped at 1 x 10^°cm~'^ and 1 x lO^^cm"^ respectively. The boundary 
conditions for the potential at the contact are set to floating, by imposing the condition that 
the potential derivatives are zero at these boundaries. The drive current vs. gate voltage is 
plotted in Fig. 12. The Si devices yield a slightly larger ballistic current compared to their 
InSb counterparts. The main reason is due to the larger density-of-states mass in Si and the 
fact that it is doubly degenerate compared to non-degenerate InSb F valley. Currently, ex- 
perimental study shows that sub 50 nm Si devices only achieve ~ 40% ballistic performance 



3l| . A lower density-of-states mass should help damp the dissipative processes to achieve 



fully ballistic performances in InSb devices and performance then could be better than Si. 



V. CONCLUSION 



Band structure of III-V material InSb thin films is calculated using empirical pseudopo- 
tential method (EPM). F valley in InSb remains the lowest lying conduction valley (a desir- 
able trait if high mobility characteristic is required) despite size quantization effects but its 
isotropic effective mass increases with decrease in the film thickness. A NEGF calculation 
of the InSb and Si double-gated devices reveals that they have comparable ballistic drive 
current. 
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Appendix A: treatment of spin-orbit coupling 



The first starting point is the formulas in Ref. 26|, such as equation (A3) in it. To 
explicitly write down the formula of SOC, we shall calculate one of the terms of Eq. (A4) 
in Ref. |26| denoted as (II) here. We consider 

{11) = -Y,{k'\bt){bt\A\k) 

t 

\^ — y^y^ I dridr2df3e-''''"'^iJnim{\ri - Rj - t^\,9i, ^i) 



t—nlm ij 



xrnim{\r2 - Rj - fd, ^2, 02)^(1^3 " R, - r.|)/(f3 - R^ - f.)(5(f2 - r^e^^''^^' (Al) 

where n, / and m are quantum numbers characterizing the core states, t = nlm, Rj is the 
point located the unit cell, are the coordinates of atoms in the basis, ^(|r|) = i^^j^i^^p^; 
with V being potential of atomic nuclei, I is the angular momentum operator, 6 is the Dirac 
6 function, Qceii is the volume of the unit cell, and Nceii is the number of unit cells in the 
crystal. By performing the integral over f^, and then let r/=ri — Rj — fi and r2=r2 — Rj — fi. 
Then we have 

t=nlm ij 

X [ I rfrl^:,„(r^, <P2)ar',)l{r^2y^-% (A2) 



We denote the first square bracket in Eq. (]A2|) by [A], the second one by [B] and the third 
square bracket by [C]. Each of these terms will now be individually analyzed. [A] includes 
the local position information of all atoms in a unit cell, as given by 

[A] = -^5^e^^^-, (A3) 
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where A k = k' — k is the change in the wave vector. When we consider single-element 
crystal, there is only one kind of atom. If we set riceii as the number of all atoms in a unit 
cell, then we have 

[A] = -^(^)E^^''-^ 

ncell i l-cell "~ 
^ ^atom 

where 

S(Ak) = ^Ye'''^'"' (A5) 
riceii^ 

is the counterpart of the structure factor in Ref. j^], and flatom is the atomic volume in a 

unit cell so that Qatom = ^cell/ncell- 

The second square bracket [B] in Eq. (IA2p is 

[B] = J2cJ drMk'r)Y;,,{e,<jy)Rni{r)Yr{eA) 

l'=0 ^ 

oo « /'27r nTV 

= Y.^i'[ r^Ji'{k'r)Rr.i{r)dr][ smeY,l,,{e,cl>)Yr{9,(^)de] (A6) 



where = (— l)''i''[47r(2/' + 1)]^''^, Y^^ are spherical harmonics and jii are spherical Bessel 
functions. We use the formula e^^'^ = X]/^o*'[^^(2^ + '^)Y^'^ii{kr)Yik and Yik{6,(f)) are the 
spherical harmonic with the rotational index m = in any coordinate system with z in the 
k direction. Because of the orthogonality of spherical harmonics, the second square bracket 
in Eq. ( ]A6I) gives Sii'6mm', i-e., SwSmo- Hence Eq. ( 1A6I) simplifies to 

[B] = i-iyn'/^B^i,,ik')Sm,o, (A7) 

where Bni,i = fi'^^^ i'[47r(2/ + l)Y^'^ji(k'r)Rni{r)r'^dr is the counterpart of B function 
defined in Ref. [26j, Rni is the radial part of core wave function. 
The third square bracket [C] in Eq. flA2p as given by 



oo „ 

= / RUr)ar)3i'ikrydr] x [(r^l^l^^'fe)] 

i'=o ^ 

oo 

= ''^^[radial part] x [angular part] (A8) 



l'=0 
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where angular part is 

(Yrm,) = {Yrmk)6u', (aq) 

Here m in spherical harmonics of angular part must be constrained by 6mo in Eq- (1A7P . The 
same core states characterized by nlm are expressed in the different coordinates. The term 
(fi|6t) in Eq. (lAip means the core states are projected into fi coordinates system in which 
z direction is in k'. Similarly in (6f|r2), they are projected into r2 (i.e. rs) coordinates, in 
which z direction is in A;. So we can write 1^™ in Eq. (lA9p as ¥1^1. Then we get the angular 
part to be j3] 

where Pi are the associated Legendre polynomials, a is the mutual angle between the vector 
k and k' . From Eq. f lATOll . fjA9ll and (ESI), we get 

where A'„j = (^-^/^ j~ ^«[4^(^2/ + l)Y/'^Rni{r)S,{r)]i{kr)r'^dr, which is similar to the "A" 
function in Ref. ( 26|). Then we can get the term (//) as 

nl ^ ' 

where the summations shall be taken over the quantum number n and /. 

The other terms in (A3) in Ref. 26| can also derive from the similar calculations. Then 
the final expression of SOC can be written as 

^fs',ks = -^^'^ ■ X k)Ktoml (A13) 

where dg's = {s'\(^\s) and 

Xaton. = {l/kk') Ti-iy[Br^l,^ik')A'^,{k)+Ar,lik')B^^^^^^ 

^-^ ^-^ ct cosa 

nl n' ' 

(AM) 

where the first term comes from the plane wave (p)-core state (c), the second term comes 
from c-p, and the third comes from c-c term. 

/>oo 

Tlnn'i= / Rn'i{r)i{r)Rni{r)r''dr (A15) 

is about the spin orbit splitting of isolated atoms. It is just determined by local wave 
functions of core states. We shall discuss two cases. 
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Case one: two different atoms in a unit cell We set 7 =A fc-r, Ai and A2 are corresponding 
to the two atoms. We get {27 1 



where A"^'"^ = (Ai ± A2)/2 are the symmetric (antisymmetric) contributions to the spin-orbit 
Hamiltonian. 

Case two: multi-type and many atoms in a unit cell If we have m species atoms in a unit 
cell, and ni, n2, ■ ■ ■ , Um are numbers of atoms for each type atom in a unit cell. We have 

Em 

m 

HiZrks = -'^^'^ ■ X ^) E ^^(^ ^)^^' (A17) 

where Aj is for j-th species, and structure factor is 5'j(A k) = n'^^ '^nj exp(z A k-Tnf), here 
summation over rij means over the atoms of j-th species. 



If we only keep the third term c-c in Eq. flA14p (as in Refs. 27H29|) and only consider the 
main contribution from the outermost / = 1 core state, i.e. p orbit. Also note the relation 

dPi/d{cosa) = 1, for / = 1, we have 

Xatom = {kk')-^Bni,i{k')Bnl,l{k)7]nnl, (A18) 

where n is selected as the quantum number characterizing the outermost p core state, and 
rinni is spin-orbit parameter for various atoms, we write it as rjatom- If we redefine the B 
function by introducing a normalization constant C as 

POO 

Bni{k) = C / jni{kr)Rni{r)rHr, (A19) 

where C is determined by the limit \im.k-^Qk^^Bni{k) = 1. And in this procedure, an 
empirical adjustable parameter /i and a describing the ratio of the spin-orbit contributions 
for free atoms A and B, see Refs. To evaluate the B function explicitly, we may 

determine the expressions of the core state wave function Rni{r) and we use the formula in 



Ref. 



32| 



'ipnlm — CnlpXplm, (A20) 



where Cnip are coefficients of expansion, and Xpim are Slater-type orbits with integer quantum 
numbers, namely 

Xpim{r, 9, 0) = Rip{r)Yr{9, 0), (A21) 
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where 

Ri^ = [(2n^p)!]-l/2^20p)"'''+^/V"'--^e-'^'-^ (A22) 
where Cnip and Qp can be found in the tables in Ref. So the functions ip in Eq. flA20p 



are 



V 

= Rip{r)Yr{e,(^), (A23) 
In terms of these functions, we can calculate the integral as 

/•oo 

Hk) = / ji{kr)Rni{r)r'^dr 
Jo 

POD 

= V'C'nip / ji{kr)Rip{r)r^dr 

= ^C„^pAr;p/f(fc), (A24) 

p 

where Nip = [(2n,p)!]-i/2(20p)"'''+^/^ and If{k) = ji{kr)r'^'p+^e-^'^''dr. The integral of 
If{k) can be easily derived when 1 = 1. It is 



K-l)!(sinKz.)- ^"'gyff+^^^^ 



where z/ = aictaii{k / dp) . If we calculate the value of the limit /c — 0, it shall be pointed out 
that the spherical Bessel functions ji{kr) can be formulated ji{x) = ^/^Ji+i/2{.x)^ where 
J;+i/2(a;) is a Bessel function of the first kind. Then we can get 



k {nip + 2)\ 
3 /-J^ip+s 

This result is used to derive the normalization constant in Eq. flA19p 



]^Jliik) = -'-^^. (A26) 

ip 
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